Non-equilibrium transport through quantum dots with 
Dzyaloshinsky— Moriya— Kondo interaction 
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We study non-equilibrium transport through a single-orbital Anderson model in a magnetic field 
with spin-dependent hopping amplitudes. In the cotunneling regime it is described by an effective 
spin-1/2 dot with a Dzyaloshinsky-Moriya-Kondo (DMK) interaction between the spin on the dot 
and the electron spins in the leads. Using a real-time renormalization group technique we show 
that at low temperatures (i) the DMK interaction is strongly renormalized, (ii) the renormalized 
magnetic field acquires a linear voltage dependence, and (iii) the differential conductance exhibits 
a voltage asymmetry which is strongly enhanced by logarithmic corrections. We propose transport 
measurements in which these signatures can be observed. 
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Introduction. — Over the past decade it has been estab- 
lished that electronic transport measurements through 
quantum dots and single molecules can be used to probe 
various coherent spin phenomena*^ One of the simplest 
models for the theoretical description of quantum dots is 
the single-orbital Anderson model, which consists of an 
energy level that can be occupied by up to two electrons. 
This energy level is coupled to electronic reservoirs with 
hopping amplitudes t a where a = L, R. In the regime 
of suppressed charge fluctuations the system can be de- 
scribed by the Kondo model, in which the effective spin S 
on the dot is coupled to the electron spins s in the leads 
via a Heisenberg exchange interaction oc S ■ s resulting 
in rich Kondo physics^ When the attached reservoirs are 
held at different chemical potentials, electrons are trans- 
ported through the Kondo dot by cotunneling processes. 
The complicated interplay between Kondo physics and 
the non-equilibrium processes has led to the development 
of various renormalization group (RG) approaches4£ 

An interesting extension of the Anderson model is ob- 
tained by allowing spin-dependent hopping amplitudes 
t" a , where <j, a' I. As it was shown by Paaske et al& 
spin-orbit interaction in materials like SiGe may lead to 
a mixing of different orbital states of the quantum dot, 
thus resulting in an effective spin dependence of the hop- 
ping amplitudes. Observable manifestations of this effect 
are voltage asymmetries of the differential conductance^ 
and the suppression of Kondo ridges at finite magnetic 
fields,! 

However, spin-orbit interaction is not the only possi- 
bility to generate a spin dependence of the hopping am- 
plitudes. Pustilnik et al& have considered quantum dots 
with an even number of electrons in a finite magnetic 
field. Tuning the Zeeman energy to the value of an or- 
bital splitting creates a pair of degenerate levels, which 
originate from different orbital states and thus naturally 
possess different hopping amplitudes. They further de- 
rived the effective Kondo model for the case of symmetric 
hoppings to the leads, which possesses an anisotropic ex- 
change interaction oc S x s x + S v s v + AS z s z . 



In this Rapid Communication we study a single-orbital 
Anderson model coupled to electronic leads via arbitrary 
spin- and lead-dependent hopping amplitudes t" a , (see 
Fig. Q}. We demonstrate that with these two prerequi- 
sites the regime of single electron occupancy of the dot 
is generally described by an effective Kondo model with 
anisotropic exchange interaction and a Dzyaloshinsky- 
Moriya-Kondo (DMK) interaction cx d - (S x s) between 
the effective spin on the dot and the electron spins in the 
leads. 

Dzyaloshinsky-Moriya (DM) interactions arise in most 
elementary form in magnetic molecules^ For example, 
dimer molecules can be modeled by two local spins cou- 
pled by a DM term, which can give rise to a deviation of 
the orientation of the magnetization from the direction 
of the external magnetic field. This situation is of par- 
ticular relevance^ for single-molecule magnets like Mn6 
and Mni2- The DM interaction results in a canting of 
the spins that is also predicted to have characteristic sig- 
natures in the transport properties^ 

In sharp contrast to these previous works, the DMK 
interaction studied here does not act within the quantum 
dot (or molecule) but rather between the spin (or pseu- 
dospin) on the dot and the electron spins in the leads. As 
the DMK interaction is antisymmetric under a left-right 
inversion of the system, it has drastic consequences for its 
non-linear transport properties, which we study using a 
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FIG. 1. (Color online) Single-orbital Anderson model in a 
magnetic field ho- The energies of the singly occupied levels 
are £d ± ho. The dot levels are coupled to the leads via spin- 
dependent hopping amplitudes t™ CT , = — t^(d ■ tr) aa i. In 
the Kondo regime of strong Coulomb repulsions U with ed = 
—U/2 the system is mapped onto the effective model Eq. |TJ. 



real-time renormalization group (RTRG) method:— We 
find that (i) the DMK interaction is strongly renormal- 
ized during the RG flow, (ii) the renormalized magnetic 
field (RMF) h acquires a component along the DMK vec- 
tor d that grows linearly with the bias voltage V, and (iii) 
the DMK interaction results in an asymmetry of the dif- 
ferential conductance G(V) ^ G(—V), which is strongly 
enhanced by logarithmic corrections. We further discuss 
how the RMF manifests itself in the correlation functions 
and the real-time dynamics of observables, and propose 
the measurement of the differential conductance in two- 
and three-terminal set-ups for its experimental detection. 

Model. — We consider a quantum dot coupled to two 
electronic reservoirs held at chemical potentials /iL,R = 
±V/2 (Fig. [1]). Electronic transport through the dot is 
mediated by two levels with level splitting ho, which are 
assumed to be well separated from all other energy levels 
of the dot. The hopping amplitudes t" a , depend on the 
spin states as well as the attached leads. In the case of 
normal Fermi-liquid leads the matrices t" a , can be chosen 
hermitian and thus parameterized as t" a , = t\5 aa i — 
t^{d a ■ a)crcr' in terms of four real numbers i^' 2 (one can 
always choose t\ > 0) and two unit vectors d a ; er denotes 
the vector of Pauli matrices. We focus on the case d a = 
d = (cf ) where i = x,y, z. 

In the Kondo regime one can eliminate the empty 
and doubly occupied states on the dot via Schrieffer- 
Wolff transformation. 6 This yields an effective Kondo- 
type Hamiltonian 
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where J a p = 8(t^ 2 l +t% l t^)/U is the coupling between 
the spin- 1/2 on the dot and the charge density n a p — 

Sfcfc'cr c l,ka c l3k'<T in the leads, and U denotes the on-site 
Coulomb repulsion. The exchange coupling between S l 
and the spin densities s l a/3 = §£ fcfcw c lka a ia> c Pk<*> in 
the leads is given by 



The couplings jH d = 8(^t^ =F tat 2 p)/U constitute the 
anisotropic Kondo model, 5 while J ( ^ 1K = 8(t^t 2 ! — 
t 2 a t\)lV is the Dzyaloshinsky-Moriya-Kondo (DMK) 
coupling between the spin on the dot and the electron 
spins in the leads ^ 

An important parameter in the model ([1]) is the angle 
between the DMK vector d and the applied magnetic field 
ho. A parallel alignment d || ho describes, for example, 
the Kondo effect in quantum dots with an even num- 
ber of electrons^ On the other hand, the perpendicular 
case d _L h is found 6 - in the effective cotunneling model 
for a Kramers doublet in quantum dots with spin-orbit 



interactions, provided the coupling between the orbital 
angular momentum and the magnetic field is neglected. 

We stress that a non-vanishing DMK term requires a 
left-right asymmetry of the hopping amplitudes, which is 
quite naturally present in nanoscale junctions, and con- 
tains the non-local reservoir spin density slr. The cor- 
responding coupling matrix J DMK is antisymmetric in 
reservoir space, i.e. J DMK oc it v . The other couplings 
J ±:d and J are symmetric and thus spanned by 1, t x , and 
t z , where t % denote the Pauli matrices in reservoir space. 
Hence in total the model (H} contains ten couplings. Al- 
though their bare values are expressed just through the 
four parameters t L ' ^, the relations between renormalized 
couplings appear to be more complex, and all ten of them 
have to be taken into account in the RG treatment. 

Scaling analysis. — The renormalization of the cou- 
plings J is governed to leading order by the poor-man's 
scaling (PMS) equations. For the model ([T]) they have 
the form (in the following all matrix operations are per- 
formed in the reservoir space) 
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where I = ln(D/A) with A being the flow parameter. 
Its initial value is given by the band width D. The 
parametrization ^fy is preserved under the RG flow gov- 
erned by Q, i.e. no new terms are generated. The PMS 
equations ([3]) also leave invariant tr J = tr Jo = 2c and 
tr[(J^) 2 - (J d ) 2 -J 2 - (J DMK ) 2 ] = 0. The latter rela- 
tion implies that when both J and J DMK are initially zero 
one recovers the isotropic Kondo model with J- 1 = J d . 
We emphasize that despite J 1 - and J DMK formally sat- 
isfy the same RG equation, they have different symmetry 
properties in the reservoir space. 

It is worthwhile to note that the magnetic field ho does 
not explicitly enter the PMS equations ([3]). Its magni- 
tude ho may appear as a cut-off scale for the RG flow, 
but its direction is irrelevant for the PMS analysis. Hence 
the PMS equations possess a higher symmetry than the 
full system (JTJ) . This allows us to arbitrarily choose the 
direction of d, which is not affected by the RG flow. The 
Eqs. ([3]) with the initial (bare) values Jo = J(A = D) 
stated after Eqs. (P) and ((2]) can be solved exactly. Here 
we focus on the most significant features of the RG flow 
and the role of the spin-charge and DMK terms. 

The parameters trJ d and trJ^, which essentially de- 
termine the flow diagram of the model ([T]), obey 



tr J d = 2c 



1 + (Jk/A) 
1 - (T K /A) 
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4c (T K /A) 
K 1 - (T K /A) 
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where the invariant Kondo temperature is defined by 
T K = D [(trJ d - 2c)/(trJ d + 2c)] 1 / 4c , and c was defined 
above. In comparison with the conventional anisotropic 
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Kondo model^ the solution ((4]) additionally depends on 
the asymmetry parameter K = -^(trJ^) 2 — 4c 2 /(tr J^), 
which is in general different from unity unless both equal- 
ities = t^ and i£ = hold. In the latter case we have 
jDMK _ q anc j j _ cons t. during the RG flow, thus we 
recover the scaling equations of Ref. H. 

In contrast, in the generic case of broken left- right sym- 
metry J DMK and J are renormalized under the RG flow 
Q and acquire non-zero values even if one of them is 
absent initially. For example, we find tr[iJ DMK r y ] oc 
V K 2 — 1 trJ^, i.e. the DMK coupling diverges for A — > 
Tk- In the following we show that precisely J DMK and J 
cause interesting observable effects. 

Renormalized magnetic field (RMF). — One of these ef- 
fects is the renormalization of the magnetic field. For its 
derivation we have to go beyond the PMS analysis. We 
use the RTRG method^^ as it is particularly suited to 
analyze non-equilibrium effects. The derived flow equa- 
tion for the RMF can be integrated with the result 




FIG. 2. Differential conductance with (solid line) and without 
(dashed line) logarithmic corrections for t\, = 0.08, = 0.04, 
fp = 0.05, 4 = °- 03 > d II h o, and h = 100 Ik- We observe 
an asymmetry under V — >■ — V, which gets strongly enhanced 
by the logarithmic corrections. Insets: The asymmetry of 
the RMF with respect to V —¥ — V can be read off from the 
resonances m dG/dV = d 2 //dV 2 . 
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h -tr[(J c - J )t z ] -d. (5) 



The subscript c denotes the exchange couplings cut off 
at the energy scale A c = max{/io,V} 3> Ik. The first 
term in (JSJ is analogous to the first-order result for the 
RMF in the anisotropic Kondo model. The second term 
originates from the renormalization of the coupling be- 
tween the dot spin and the electron density in the leads, 
which is generated by the DMK interaction [see Eq. fl3"])]. 
Hence the linear voltage dependence of ([5]) is a conse- 
quence of the left-right asymmetry and spin dependence 
of the hopping amplitudes. We stress that the voltage 
dependence already appears in first order in the renor- 
malized exchange couplings J c , and that the effect on the 
magnitude h of the RMF is maximal in the parallel case 
d || Hq. The voltage dependence implies that a RMF will 
be generated even in the absence of an external field ho 
and that h(V) ^ h(—V). This is in sharp contrast to the 
anisotropic Kondo model^, where the voltage dependence 
only enters through logarithmic corrections in 0(J 2 ) and 
the RMF remains symmetric under inversion of the bias. 

The result ([5]) has an immediate consequence for the 
dynamics^ (S x (t) ± iS y (t)} oc e ±lht of the transverse 
magnetization. In the parallel case the precession fre- 
quency h depends on the relative sign of ho and V, which 
can be switched externally. The comparison of the two 
values of h directly yields the voltage-dependent term in 
([3]). Furthermore, the dynamical spin susceptibility pos- 
sesses^ characteristic features at = h,h ± V, which 
can also be used to probe the RMF. The experimental 
observation of the precession frequency as well as the spin 
susceptibility does, however, require time-resolved mea- 
surements on the quantum dot. Although time-resolved 
experiments have been previously performed, 17 it would 
be rather challenging to extract quantitative results for 
the RMF. Thus we propose two alternative ways to ob- 
serve the effect of the RMF in stationary quantities. 



Differential conductance. — As the DMK interaction 
breaks the inversion symmetry of the system, it has dras- 
tic consequences on its non-equilibrium transport prop- 
erties. In Fig. [2] we plot the differential conductance 
G(V) = dl/dV through the dot for a parallel alignment 
of d and ho- For |V| < h the conductance is dominated 
by elastic cotunneling processes. At V f=s ±/i the con- 
ductance possesses jumps due to the onset of inelastic 
cotunneling processes. The magnitude of these jumps 
is asymmetric with respect to V — > — V already in or- 
der J 2 (dashed line in Fig. ^ , which qualitatively agrees 
with previous perturbative results. 6 In leading order this 
asymmetry is proportional to tr[iJ,P MK J c L r z ](d ■ ho)/ho- 
In addition, the RTRG analysis shows that the zero- 
temperature broadening of these jumps is given by the 
transverse spin relaxation rates T2{V = h) and T2(V = 
—h), respectively, which differ from each other by cx 
tr[iJ c DMK J c L r 2; ](d ■ h ). This explicitly shows that the 
asymmetry of G(V) is a direct consequence of the DMK 
interaction and it can be quantified in terms of the corre- 
sponding renormalized coupling J,P MK . The asymmetry 
further depends on the angle between d and ho, taking 
its maximal value in the parallel case and vanishing in 
the perpendicular case. 

A more thorough RTRG analysis in the parallel case 
including the leading logarithmic corrections in 0(J<?) 
reveals that the magnitudes of the jumps at V = ±h ac- 
quire logarithmic terms oc ln(A c /|V =F h + ir2(T^ =F h)\) 
broadened by the corresponding transverse spin relax- 
ation rates. Due to the difference of the respective pref- 
actors the asymmetry of G(V) acquires a strong enhance- 
ment (see the solid line in Fig. [5]), which even influences 
the elastic cotunneling regime \V\ < h. 

The logarithmic resonances at the onset of inelastic 
cotunneling at V = ±h correspond to the extrema of 
the derivative of the differential conductance dG/dV — 
d 2 I/dV 2 . Thus the RMF h(^fV) can be read off from 
the positions of its dip and peak respectively (see insets 
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FIG. 3. (Color online) a) Sketch of the three-terminal set-up. 
The stationary current /(V) flows between the left and right 
leads while the probe lead is used to measure the RMF. b) 
Conductance G p in the probe lead for J p = 10 -3 and ho = 
2V. All other parameters are as in Fig. [2] c) The RMF © 
can be read off from the resonances in d/dy [G p (/i p + V/2)]. 
The curves shown correspond (from top to bottom) to V = 
20 T K , 40 T K , . . . , 160 T K . The dashed line is a guide to the 
eye indicating the linear voltage dependence of h(V). 

in Fig. [5]). The second term in (J5J) is given by the sum 
of the two positions. We note that due to the restriction 
of our quantitative analysis to the weak-coupling regime 
the asymmetry h(V) ^ h(—V) observed in Fig.[3]is rela- 
tively small. However, as it turns out to be a generic con- 
sequence of asymmetric and spin-dependent hopping am- 
plitudes we expect it to be present in the strong-coupling 
regime as well. 

Three-terminal set-up. — Another way to experimen- 
tally observe the linear voltage dependence of the RMF is 
provided by the use of the three-terminal set-up sketched 
in Fig. [3]a. A similar set-up has been used by Leturcq et 
al*£ to measure the splitting of the Kondo resonance in 
the spectral density. Specifically, the three-terminal set- 
up consists of the system (p} with an additional, weakly 



coupled probe lead modeled by J2k<r( e k ~ Mp) c lk* c pkv 
with exchange interaction J p J2 a =L r & ' ( s ap + s pa) + 
J p S ■ s pp , where |J P | -C \J%p\, \J a p\- We assume sym- 
metric and isotropic exchange couplings to the probe lead 
and note in passing that the calculation can easily be gen- 
eralized to more complicated situations. As the probe 
lead is only coupled weakly, its influence on the PMS 
equations and the RMF is negligible, and the conduc- 
tance G p = dlp/dfip can be calculated^ using perturba- 
tion theory in J p . The result is shown in Fig.[3]b: G p pos- 
sesses characteristic steps at /z p = ±(h— V/2), ±(h+V/2) 
which are broadened by the transverse spin relaxation 
rate. When plotting dG p /dV as a function of fj, p + V/2 
one can directly read off the RMF from the position of 
the first resonance. The results shown in Fig. [3jc clearly 
reveal the linear voltage dependence of the resonance po- 
sition and thus the RMF. 

Conclusion. — We have studied a Kondo dot with ad- 
ditional DMK interaction, which is amplified during the 
RG flow. The DMK term results from the broken in- 
version symmetry of the system and manifests itself in 
the non-equilibrium transport properties. In particular, 
we showed that the RMF acquires a component along 
the DMK vector depending linearly on the bias voltage, 
which can be experimentally detected in two- and three- 
terminal set-ups. In addition, the differential conduc- 
tance becomes asymmetric with respect to the bias. As 
these fingerprints of the DMK interaction rely on broken 
spin and spatial symmetry, they are also expected be- 
yond the weak-coupling regime. A complete description 
of the strong-coupling regime remains, however, an open 
problem. 
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